Mon. Not. R. Astron. Soc. 000, [TV?? (2002) Printed 11 July 2008 (MN WT^i style file v2.2) 



A Detailed Study of 2S 0114+650 with the RXTE 



S. A. Farrell 1 ' 2 *, R. K. Sood 2 , P. M. O'Neill 3 ' 4 and S. Dieters 5 

1 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, 9 Avenue du Colonel Roche, 31028 Toulouse Cedex 04, France 
^School of F 'EMS, UNSW&ADFA, Canberra ACT 2600, Australia 

3 Astrophysics Group, Imperial College, Prince Consort Road, London SW7 2AZ, UK 

^School of Computing & Mathematics, Charles Sturt University, Locked Bag 588, Wagga Wagga NSW 2678, Australia 
School of Mathematics and Physics, University of Tasmania, Private Bag 21, Hobart Tasmania 1001, Australia 



Released 2007 Xxxxx XX 



ABSTRACT 

We present the results of a detailed study of the high mass X-ray binary 2S 0114+650 
made with the pointed instruments onboard the Rossi X-ray Timing Explorer. The 
spectral and temporal behaviour of this source was examined over the pulse, orbital, 
and super-orbital timescales, covering ~2 cycles of the 30.7 d super-orbital modulation. 
Marginal evidence for variability of the power law photon index over the pulse period 
was identified, similar to that observed from other X-ray pulsars. If this variability is 
real it can be attributed to a varying viewing geometry of the accretion region with 
the spin of the neutron star. Variability of the neutral hydrogen column density over 
the orbital period was observed, which we attribute to the line of sight motion of the 
neutron star through the dense circumstellar environment. A reduction in the power 
law photon index was observed during the orbital maximum, which we speculate is 
due to absorption effects as the neutron star passes behind a heavily absorbing region 
near the base of the supergiant companion's wind. No significant variability of the 
column density was observed over the super-orbital period, indicating that variable 
obscuration by a precessing warp in an accretion disc is not the mechanism behind the 
super-orbital modulation. In contrast, a significant increase in the power law photon 
index was observed during the super-orbital minimum. We conclude that the observed 
super-orbital modulation is tied to variability in the mass accretion rate due to some 
as yet unidentified mechanism. 

Key words: accretion, accretion discs - stars: neutron - X-rays: binaries - pulsars: 
individual: 2S 0114+650 



1 INTRODUCTION 

The high mass X-ray binary 2S 0114+650 has been 
shown to be a rather unusual system. It exhibits prop- 
erties consistent with both Be and super-R iant X-ray bi- 
naries (|Crampton. Hutchings fc Cowlevlll985l ). with signif- 
icant broadband temporal and spectral variability over a 
wide range of timesc ales (e.g. iKoenigsberger et al.l 1 19831; 
ICrampton et al.1 Il985l h The derived low X-ray luminosity 
(L a ~1.1 x 10 36 er g s" 1 in the 3 - 20 keV band for a 
distance of 7.2 kpc; lHall et ail |2000T ) implies that spheri- 
cal a ccretion takes place via th e stellar wind of the donor 
star (|Li fc van den Huevellll999h . However, the Hell 4686 A 
line (a common signature of the presence of an accretion 
disc) ha s been observed weakly in emiss i on on a few oc- 
casions l|Aab. Bvchkova fc Kopylovl Il983l : ICrampton et al.1 
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Il985l : Ivan Kerkwiik fc Water j|l989h . possibly linked to the 
presence of a small transient accretion disc. 

Coherent X-ray (and possibly faint optical) pulsations 
at a period of ~2.7 h h ave been shown to be a persis- 
tent feature of the system (|Finlev. Belloni fc Cassinellilll992l : 
iTavlor et ai]|l995t) and have been variously attributed to (3 
Cephei-type pulsations/tidal oscill ations of the dono r star 
and the neutron star spin period. ISood et al.1 (|2006T ) anal- 
ysed the evolution of the pulse period in ~9.5 yr of Rossi 
X-ray Timing Explorer (RXTE) All Sky Monitor (ASM) 
data, identifying two episodes of torque reversal. Their con- 
clusions support the scenario that 2S 0114+650 contains a 
super-slow rotator - one of the slowest spinning X-ray pul- 
sars yet discovered. 

Vari ability at ~11.6 d ha s been observed in both 
optical dCrampton et alj 19851 ) and X-ray wavelengths 
l|Corbet. Finlev fc Peeld I1999TI. consistent with the orbital 
period of the system. Crampton et al.l l|l985h were the first 
to confirm the binary nature of 2S 0114+650 using radial 
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velocity measurements (with a period of 11.588 ± 0.003 d 
for an eccentricity of 0.16 ± 0.07, and 11.591 ± 0.003 d for 
an eccentricity of 0.0). Although they could not distinguish 
between a circular and an eccentric orbit, they concluded 
that the orbital motion was most likely eccentric as the sys- 
tem had presumably gone throug h a supernova stage which 
should lead to a perturbed orbit. ICorbet et alj (|l999h anal- 
ysed ~2.5 yr of RXTE ASM observations of 2S 0114+650 
between 1996 - 1999, identifying X-ray modulations close 
to both the optically derived orbital period and the pro- 
posed pulse period. Fitting a sine wave to the ASM light 
curve yielded a period of 11.63 ± 0.007 d. The inconsistency 
between the optical and X-ray orbital period values was ex- 
plained as either due to variable contamination of absorption 
lines by low equivalent width emission lines (resulting from 
the X-ray irradiation of the donor star) or possibly variabil- 
ity in the X-ray orbital modulation. 

The determination of the orbital period was revis- 
ited recently in b o th op tical and X-ray wavele ngths by 
IGrundstrom et al.l (120071) and |Wen et al.l (|2006h respec- 



tively. IGrundstrom et al] l|2007 ) performed a long-term 
spectroscopic monitoring program of the Ha emission in 2S 
0114+650 and presented revised orbital elements with a pe- 
riod of 11.59 83 ± 0.000 6 d an d a non-zero eccentricity of 
0.18 ± 0.05. IWen et al.l ([20061 ) re-analysed the ASM data 
and reported a new value for the orbital period of 11.599 ± 
0.005 d, consistent with the new optical period. 

A super-orbital modulation at 30.7 ± 0.1 d was reported 
by IFarrell. Sood fc O'Neilll (|2006|) from the analysis of 8.5 
yr of ASM data. IWen et al.l (|2006l ) subsequently reported a 
new impro ved value of 3 0.75 ± 0.03 d from re-analysis of the 
same data. ISood et"aL I 1120061) analysed the evolution of the 
super-orbital period, finding it to be stable and concluding 
that 2S 0114+650 is only the fourth reported system to show 
stable super-orbital variability. 

The standard model invoked to explain super-orbital 
variability in X-ray binaries is periodic obscuration or re- 
flection via the precess ion of a warped accretion disc (e.g. 
lOgilvie fc~ Dubus 200ll). The apparent lack of a persistent 
disc in this system thus raises doubts about the validity of 
this model as an e xplanation for the observed 30.7 d mod- 
ulation. Recently, iKoenigsberger et al.l l|2006l ) presented a 
model by which both the pulse and super-orbital modula- 
tions in this system could be explained by tidally induced os- 
cillations in the donor star. However, the precise mechanisms 
by which these oscillations may translate into a structured 
stellar wind are not clear. In addition, super-orbital variabil- 
ity is predicted for circular orbits while the eccentric orbit 
models predict strong variation s on orbital timescales asso - 
ciated with periastron passage (|Koenigsberger et al.| [2006). 
These results are thus inc onsistent with the non-zer o eccen- 
tricity recently derived by IGrundstrom et al] (|2007l ). 

In this paper we present the results of a series of X- 
ray observations of 2S 0114+650 perf ormed with the Pro - 
portional Counter Array (PCA; see IJahoda et al] Il996l ) 
and High Energy X -ray Timing Experiment (HEXTE; see 
iGruber et al.lll996n instruments onboard the RXTE satel- 
lit e. The prelim i nary results of these studies were published 
m IFarrell et al] (2007). The temporal and spectral variabil- 
ity of this source was studied in detail over the pulse, orbital 
and super-orbital timescales. In § 2 the data reduction steps 
that were taken are explained. A description of the analysis 



methods employed and the results obtained are provided in 
§ 3. A discussion of our results is provided in § 4, and § 5 
presents the conclusions that we have drawn. 



2 DATA REDUCTION 

Between 2005 May 15 - 2005 June 14 (MJD 53505 - 53535) 
and 2005 December 13 - 2006 January 12 (MJD 53717 - 
53747) 2S 0114+650 was observed for ~220 ks with the PCA 
and HEXTE instruments. Two separate runs of 11 observa- 
tions were made, with each individual observation ~3 hr 
in duration so as to cover one complete cycle of the pulse 
period per observation. Each run covered one super-orbital 
cycle. The two runs were spaced approximately six months 
apart, with observations scheduled every few days. Table [1] 
lists the 22 observations with the date and exposure for each 
observation. 

The Standard2f mode PCA data (16 s integration time, 
128 energy channels) were used for our timing and spectral 
analyses, as time-resolution less than 16 s was not required. 
The HEXTE Archive mode data (16 s integration time, 64 
energy channels) were used for our spectral analyses to in- 
vestigate the high energy end of the spectrum. The FTOOLS 
v6.1 software package was used for the reduction and anal- 
ysis steps outlined below. 

Filtering criteria were employed to exclude times when 
the pointing direction was less than 10° from the Earth 
horizon, the offset from the target direction was greater 
than 0.02° or the electron contamination fraction exceeded 
0.1. Data obtained within 30 min of the spacecraft passage 
through the South Atlantic Anomaly (SAA) were also ex- 
cluded to remove the effect of residual high-background lev- 
els due to induced radioactivity. 

Spectra were extracted from all anode layers for each 
Proportional Counter Unit (PCU) in the PCA from all data 
sets, excluding PCU0 due to the increase in the background 
that occurred therein with the loss of the propane layeiQ. 
Spectra were also extracted for both HEXTE clusters (clus- 
ter and 1) from all data sets. Due to the intermittent 
failure of the rocking mechanism in cluster during our sec- 
ond set of observation^, background data for this cluster 
are not available for observations 12 - 19. As a result, clus- 
ter data from these observations are not included in our 
spectral analyses. 

Three light curves covering the 2-9 keV (A), 9-20 
keV (B) and 20 - 40 keV (C) energy ranges were extracted 
from all layers of PCU1, PCU2, PCU3 and PCU4 separately. 
The energy boundaries were chosen in line with those of the 
standard product light curves provided with the data. As 
the average count rate before background subtraction was 
less than 40 counts s _1 , the 2005 November 28 mission- long 
faint source model was used to generate PCA background 
spectra and light curves. HEXTE background spectra for 
each cluster were extracted from the "off-source" data, when 
the cluster was aimed in a direction offset from the pointing 
direction. PCA response matrices were generated for each 



1 http:/ /universe. gsfc.nasa.gov/xrays/programs/rxte/pca/doc/ 
bkg/bkg-2002/ 

2 http:/ /heasarc. gsfc.nasa.gov/docs/xte/whatsnew/big. html 



A Detailed Study of 2S 0114+650 with the RXTE 3 



Table 1. RXTE pointed observations of 2S 0114+650. 



Obs. 


Date 
(MJD) 


Exposure 
(s) 


Obs. 


Date 
(MJD) 


Exposure 
« 


1 


53506 


7,104 


12 


53717 


7,216 


2 


53509 


7,488 


13 


53721 


10,176 


3 


53513 


5,856 


11 


53723 


8,608 


4 


53515 


7,408 


15 


53726 


8,064 


5 


53517 


9,120 


16 


53729 


8,336 


6 


53520 


9,616 


17 


53732 


1,632 


7 


53523 


8,960 


18 


53735 


7,888 


8 


53526 


6,640 


19 


53738 


7,616 


9 


53529 


7,152 


20 


53741 


7,664 


10 


53533 


7,056 


21 


53745 


6,752 


11 


53535 


7,408 


22 


53747 


8,592 



Table 2. XSPEC spectral models applied to the combined PCA 
and HEXTE spectra. A photo-electric absorption component 
(WABS in xspec) was included for all models. 



Model 


XSPEC Models 


X 2 /DOF 


Power law + HE cut-off 


POW*HIGHE 


28.8/80 


Black body + bremsstrahlung 


BBODY+BREMS 


29.6/80 


Cut-off power law 


cut-off 


32.5/81 


Power law + bremsstrahlung 


POW+BREMS 


33.2/80 


Comptonisation 


COMPTT 


56.9/79 


Thermal bremsstrahlung 


BREMS 


63.1/82 


Comptonisation 


COMPST 


94.7/81 


Black body 


BBODY 


3890/82 



individual source spectrum. The 2000 May 26 HEXTE re- 
sponse and ancillary response matrices were used, as the 
response of each cluster has not changed since that time. As 
the observed count rates were relatively low, a deadtime cor- 
rection was not applied to either the spectra or light curves. 

The individual PCU spectra were then combined 
(source and background, separately) to maximise the signal- 
to-noise ratio. The PCA response matrices were also com- 
bined, weighted by the exposure values for each PCU__. For 
the analysis of the complete spectrum covering all obser- 
vations, a 1% systematic errou was added to the resulting 
spectrum to account for uncertainties in the response ma- 
trix such as the feature present at ~5 keV due to the xenon 
L-edges. For the phase resolved spectral analyses no sys- 
tematic error was applied as our aim was to constrain the 
spectral parameters to the highest level possible, particu- 
larly in the low energy range of the spectrum. Instead, the 
features around ~5 keV were modelled using a multiplica- 
tive absorption edge component during the spectral fitting. 

At energies below 3 and above 30 keV there were poor 
counting statistics and so these data were excluded from the 
PCA spectrum. The observed scatter in the HEXTE spectra 
is of the lev el expected for systema tic uncertainties in the 
background (|Rothschild et al1ll998l ). and so a 1% system- 
atic uncertainty was added to its PHA channel data. Ener- 
gies less than the 17 keV low energy threshold were also ex- 
cluded in the HEXTE spectra. The low flux of 2S 0114+650 
at high energies resulted in poor statistics above 50 keV in 
the HEXTE spectra. As such, energies above this range were 
excluded for our analyses. The individual PCU light curves 
were averaged together to give source and background light 
curves with the best signal-to-noise ratio. Barycentric cor- 
rection was applied and background subtraction performed. 



3 DATA ANALYSIS 

After the data reduction steps outlined above, spectral 
models were fitted to the combined PCA and HEXTE 
(cluster and 1, separately) spectra in the 3-50 keV energy 
range using xspec vl2.3. Eight models commonly applied 



3 http: / /heasarc. gsfc.nasa.gov/docs/xte/recipes/pcu_combine. html 

4 http: / /heasarc . gsfc . nasa.gov /docs/ xte / ftools / xtefaq_answers .html 



Table 3. Best fit X-ray spectral parameters of 2S 0114+650, 
achieved using an absorbed power law with a high energy expo- 
nential cut-off. The flux quoted is not corrected for absorption. 



Spectral Parameter 


Value 


Units 


Column Density 


(3.2 t°o 9 s ) x 10 22 


atoms cm -2 


Photon Index 


1.1 ± 0.1 




Cut-off Energy 


6.0 ± 0.7 


keV 


Folding Energy 


15 tl 


keV 


Flux (3 - 50 keV) 


2.3 x lO" 10 


ergs cm" 2 s~ 1 



to X-ray binaries were tested, with the resulting x 2 values 
and degrees of freedom presented in Table [2] The best fit 
model was found to be an absorbed power law with a high 
energy (HE) exponential cut-off, with parameter values 
consistent with the previou s ly ob t ained values listed in the 
litera t ure (Yamauc hi et al.l 1990; Apparao, Bis ht fc Singhl 
19911; iFinlev. Taylor, fc Belkmil Il99l; lHall et all 
Koenigsberger et al.l 20031; Bonning fc Falanga 



20061 : iMukheriee fc Paull l2006l 



Kocnigsbcrgcr et al.l I2UUJI; IfcSonnmg fe falangal zwb; 
Filippova et all 120051: Iden Hartog et aubood: iMasetti et all 



The reduced x values obtained by fitting 5 of the 
7 remaining models were < 1, also indicating acceptable 
fits. Many of the models (such as the CUTOFF, COMPTT 
and COMPST) have the same general shape as the empir- 
ical power law + HE cut-off model, but attempt to pro- 
vide a better physical description of the underlying emission 
mechanisms. The other models (the black body and ther- 
mal bremsstrahlung) provided poor fits on their own, and 
required more complicated combinations in order to ade- 
quately model the overall spectrum. For these more com- 
plicated models we applied the F-test in xspec in order to 
justify the addition of the extra components. The addition of 
a soft black body component to the thermal bremsstrahlung 
model produced an F-statistic of 12.9 and a probability of 5 
x 10~ n , indicating a significant improvement in the fit. The 
addition of a power law componen t to the bremsstrahlung 
model (a combination suggested by Iden Hartog et al.ll2006f i 
yielded an F-statistic of 10.3 and a probability of 4 x 10 -9 , 
also indicating an improved fit. 

While we cannot discriminate between 6 of the 
models from spectral fitting alone, we have adopted 
the absorbed power law with a HE cut-off model 
(WABS*POWER*HIGHE) as it is commonly used to em- 
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Figure 2. PCA light curves of observations 1-11 (left) and 12 - 22 (right): 2-9 keV (bottom panels), 9-20 keV (middle panels) and 
20 - 40 keV (top panels). 
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Figure 3. Best fit spectral 



(blue diamonds; lYamauchi et al.l 


19901; lApparao et al.l 


1991; 


Finlev et. all Il994l: 


Hall et al. 2000 


: Kocniesbereer et al.l 


2003: 


Bonning & Falanga 


20051; iFilippova et alj|2005l;lden Hartoe et al. 


20061: iMasetti et al. 


2008: iMukheriee & Paulll200dl. From toD to 



bottom: neutral hydrogen column density, power law photon in- 
dex, cut-off energy and folding energy. 



episodes can be seen during many of the observations oc- 
curring in phase with the maximum of the pulse period 
(see e.g. the light curve from observation 13 in Figure 
IA3I in the Appendix) . These flaring episodes appear to be 
pseudo-periodic, possibly indicati ve of the 850 s , 894 s , 
or ~2,000 s periods repor ted by lYamauchi et al.l lll990h. 
iKoenigsberger et alj (|l983l ) and lApparao. Bisht fc Singhl 
1 199ll ) resp ectively. In addition to these reports of periodic 
modulation. IMasetti et all (|2006t ) reported observing "X-ray 
flickering" on timescales of minutes to hours, which they 
speculated were a result of random inhomogeneities in the 
accretion flow onto the neutron star. 

Two distinct flares were seen to occur during observa- 
tion 6 (Figure I Al|) despite the otherwise low flux from the 
source at this time. These flares are predominantly present 
at low energies, indicating a source within the PCA field 
of view a s the origin and ruli ng out gamma ray bursts as 
the cause. IMasetti et al.l (2006) suggested that similar flares 
observed in the BeppoSAX data were due to anisotropies 
in the density of the stellar wind (and thus the accreting 
material). 

The light curve of observation 22 (Figure [A4|l exhibits 
an episode of very low flux during the first hour of the obser- 
vation, after which the pulse period materialises with large 
amplitude. As mentioned earlier, all data where the pointing 
offset was greater than 0.02° from the target source were 
discarded. As such this low-flux episode is not due to the 
spacecraft slewing into position. 

Spectral and timing analyses were performed over the 
pulse, orbital and super-orbital period timescales. Analyses 
of the flaring episodes observed at the peak of the pulse 
period were also performed. These analyses are described 
below. 



pirically describe the shape of the spectra of accreting X-ray 
pulsars, and has b een used in prev ious spectral analyses of 
2S 0114+650 (e.g. lHall et aLll2000l ). 

Inclusion of an Fe Ka Gaussian emission line at 6.4 keV 
did not improve the fit significantly in the overall spectrum. 
Similarly, no evidence was found of the ~22 keV cyclotron 
resonance scattering fe ature (CRSF) tenta t ively identified 
in INTEGRAL data bv lBonning fc Falangal d2005h. Ou r re- 
sults thus suppor t the fi ndings of IMasetti et al.l " i[2006h and 
Iden Hartog et all |2006), who failed to detect any evidence 
of cyclotron features in the spectrum of 2S 0114+650 us- 
ing data from the BeppoSAX and INTEGRAL satellites re- 
spectively. Upper limits at the 95% confidence level for the 
equivalent width of the iron line and the depth of the 22 
keV cyclotron feature of 55 eV and 0.3 respectively were 
obtained. 

Figure [T] shows the combined PCA and HEXTE spectra 
from all observations fitted with the best fit model described 
above. Figure [5] shows the PCA A, B and C light curves 
from both observation runs covering MJD 53505 - 53535 
and MJD 53717 - 53747. Figure [3] compares the best fit 
spectral parameters with values from previous observations 
listed in the literature. The PCA A band light curves for all 
22 individual observations are included in Appendix A. 

The ~2.7 h pulse period can clearly be seen in most 
of the light curves, appearing as a smooth near-sinusoidal 
modulation (see Figures IA1I - IA4[I . Intermittent flaring 



3.1 Pulse Variability 

An epoch folding \ 2 search was performed on the PCA 2 - 
9 keV light curve using the ftools task efsearch in the 
9,360 - 9,720 s range with 0.36 s resolution. Light curves 
from the two observation runs (observations 1-11 and 12 - 
22) were analysed separately to determine whether the pulse 
period varied significantly over the 6 month period separat- 
ing the observation runs. The values for the pulse period 
were found to be 2.65 ± 0.01 h and 2.644 ± 0.003 h respec- 
tively, indicating the pulse period did not vary significantly 
over the 6 month period between the observation runs. The 
efsearch analysis of the light curve covering both runs of 
observations (with a reduced resolution so as to "smooth 
out" aliasing effects) yielded a best fit period of 2.650 ± 
0.003 h FWHM. Figures H [5] and |6] show the \ 2 vs period 
plots from the efsearch analysis. 

In order to investigate the profile of the pulse period 
the A, B and C (Figure [2]) PCA light curves covering all 
observations were folded over the best fit 2.65 h period us- 
ing the task efold, with phase set at the flux minimum 
at MJD 53504.912 and with the errors being the standard 
error for each phase bin. The pulse profile was found to be 
approximately sinusoidal with a short "spike" around phase 
0.2, a larger feature appearing around phase 0.85, and little 
variation between energy bands (Figure [7} . 

If 2S 0114+650 contains a super-slow X-ray pulsar with 
a 2.65 h spin period, the pulse period should exhibit Doppler 
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Figure 4. Epoch folding \ 2 search for the pulse period in the 2 
- 9 keV PCA light curve covering the first observation run (MJD 
53505 - 53535). 
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Figure 6. Epoch folding x 2 search for the pulse period in the 2 
— 9 keV PCA light curve covering all observations (MJD 53505 - 
53747). 
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Figure 5. Epoch folding \ 2 search for the pulse period in the 
2 — 9 keV PCA light curve covering the second observation run 
(MJD 53717 - 53747). 



shifting over the 11.6 d orbital period by ~16 s. Likewise, if 
the 30.7 d super-orbital period is related to the orbital period 
of a second compa ct object in the system (as postulated by 
iFarrell et al.ll2006l ), the pulse period should exhibit Doppler 
shifting over this timescale by ~12 s (the orbital angle of 
inclination permitting). As each of the 22 observations cover 
less than two pulse cycles, we were unable to measure the 
pulse arrival times in each data set to the required level 
of accuracy to constrain any variability over the orbital or 
super-orbital timescales. Additionally, it is unclear whether 
it is possible to measure the pulse period to the required 
level of accuracy (~0.1% of the pulse period) with further 
observations. 

To investigate the spectral variability over the pulse pe- 
riod, we extracted ten phase resolved spectra separated by 
0.1 phase bins over the 2.65 h pulse period using the same 
method outlined earlier. For these analyses, only the PCA 
3-30 keV data were utilised due to the poor statistics in 
the HEXTE data resulting from the relatively short expo- 
sures in each phase bin. As the objective of these analyses 
was to constrain the spectral parameters in each phase bin 
as tightly as possible, the 1% systematic error that was ap- 
plied to the combined spectrum from all data was removed. 
The incorrect modelling of the xenon L-edge around ~5 keV 




Pulse Phase 

Figure 7. PCA Standard2 A (top), B (middle) and C (bottom) 
light curves of all observations folded over the 2.65 h pulse period, 
with error bars determined from the standard error. Two cycles 
are shown for clarity. 
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Figure 8. The effects of varying spectral parameters on the ratio 
between the best fit model and the phase resolved spectra from 
simulations: a) varying the column density by ± 1 X 10 22 cm -2 , 
b) varying the photon index by ± 0.1, c) varying the cut-off energy 
by ± 1 keV, and d) varying the folding energy by ± 3 keV. 



was instead accounted for by an additional absorption com- 
ponent for these analyses, as outlined in §2. 

Comparing the individual phase resolved spectra with 
the best fit model obtained using the overall spectrum in 
Figure [T] allows for the identification of any energy specific 
spectral variability. In order to determine the source of any 
spectral variability, we first performed simulations using the 
overall spectrum. The spectral parameters of the best fit 
model were modified and the ratio between the modified 
model and the overall spectrum were calculated. As can be 
seen in Figure [8] varying each of the spectral parameters 
in turn resulted in quite specific patterns within the ratio 
plots. Varying the power law normalisation parameter (not 
shown in Figure [8} simply shifted the ratio values up or 
down, while keeping the slope of the ratios flat. Varying the 
column density modified the ratios only below ~5 keV, while 
varying the parameters of the high energy cut-off (i.e. the 
cut-off energy and folding energy parameters) modified the 
ratios above ~5 keV. Varying the photon index modified the 
slope of the ratio lines, while varying all parameters simul- 
taneously produced a combination of all above mentioned 
effects. 

The ratios between each of the pulse phase resolved 
spectra and the best-fit model obtained from the overall 
spectrum were then plotted in order to identify any spec- 
tral variability over the pulse period, and to identify the 
energy ranges in which the spectra were variable. The re- 
sults are plotted in Figure [5] For clarity the spectral model 
was normalised for each of the spectra through the addition 
of a constant parameter, removing the offset produced by 
the differing power law normalisation. 

The ratio plots shown in Figure [9] indicate significant 
spectral variability over the pulse period. Comparison with 
the results of the simulations presented in Figure[8]indicates 
likely variability in the slope of the power law and in the 
column density parameters, with no evidence for variability 
in the high energy cut-off parameters. Additional features in 



the ratios around ~6 keV indicate possible flux dependant 
variability in the Fe Ka emission. 

Based on this comparison, we fit each of the pulse phase 
resolved spectra simultaneously with an absorbed power law 
model with a high energy cut-off, with the cut-off energy 
and folding energy parameters tied together, but with the 
column density and power law parameters free to vary in- 
dependently. In this method the cut-off parameters are not 
fixed, but are linked so that while they were free to vary 
overall, the values were the same for each of the phase re- 
solved spectra. 

Significant features in the residuals were observed 
around ~5 keV and ~6 keV, which we attribute to the 
incomplete modelling of the xenon L-edges and to fluores- 
cent Fe Ka emission. Additional model components com- 
prising an absorption edge around ~5 keV and a Gaus- 
sian emission line at 6.4 keV (with the energy frozen at 
6.4 keV but the sigma and normalization parameters free 
to vary between the spectra) were thus included. While we 
cannot determine whether the Fe line is a result of flu- 
orescent emission from 2S 0114+650 or from the diffuse 
galactic background, similar features have been seen in the 
spect r um before (e.g . jYamauchi et al] Il990l; I Apparao ct al. 
19911; lEbisawal Il997l; iHall et al.l 120001 ; iMasetti et all l2006h 
Mukheriee fe Paull 120061 ) and so we assume that the ob- 
served line is related to the source. The reduced \ 2 statistic 
obtained from this fitting process was 0.531 with 478 degrees 
of freedom, indicating a very good fit. 

The column density and power law photon index pa- 
rameters were then determined and plotted against pulse 
phase and compared with the PCA 2-9 keV light curve 
folded over the same period (Figure [10)) . The reduced \ 2 
value for the hypothesis of parameter constancy was then 
calculated, with values of 1.41 (null hypothesis probability of 
18%) and 1.23 (null hypothesis probability of 27%) obtained 
for the column density and photon index respectively, pro- 
viding marginal evidence of variability in both parameters 
over the pulse period. The linear correlation coefficient was 
also calculated for each parameter compared to the folded 
pulse profiles. The value obtained for the column density was 
-0.22, indicating a weak anti-correlation, while the value ob- 
tained for the photon index was -0.82, indicating a strong 
anti-correlation. While the addition of an Fe Ka emission 
line was required for the best fit, we were unable to con- 
strain the equivalent width of the line sufficiently in order 
to determine whether it is variable over the pulse period. 
The pulse phase resolved spectra are shown in Figures IB1I 
and lB2l in the Appendix. 

3.2 Orbital Variability 

Approximately 6 orbital cycles were sampled during the 
course of our RXTE observations. Due to the limited num- 
ber of cycles and relatively poor phase coverage, the value of 
the orbital period could not be accurate ly constrained. As 
such, the value of 11.599 d obtained by I Wen et al] (|2006f) 
was used for our analyses. The A, B and C (Figure ITT]) PCA 
light curves covering all observations were folded over 11.599 
d, with phase set at the flux minimum at MJD 53503.7. 
The profile of the modulation does not appear to vary sig- 
nificantly with energy. 

Using the same method outlined above for the pulse 
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Figure 9. Ratios between the pulse phase resolved spectra and the best fit model from the overall spectrum (top inset figures) plotted 
against the folded pulse profile (histogram). For clarity only representative error bars are plotted on the ratio plots. 



period, we extracted five phase resolved spectra separated 
by 0.2 phase bins from the PCA data over the 11.599 d 
orbital period. The ratios between each of the orbital phase 
resolved spectra and the best-fit model obtained from the 
overall spectrum were then plotted in order to identify any 
spectral variability. The results are plotted in Figure [12] 

The ratio plots shown in Figure [12] indicate significant 
spectral variability over the orbital period similar to that 
seen over the pulse period. Comparison with the results of 
the simulations presented in Figure [8] indicates likely vari- 
ability in the slope of the power law and in the column den- 
sity parameters, with no evidence for variability in the high 
energy cut-off parameters. Features in the ratios around ~6 
keV are present similar to that seen in the pulse ratio plots, 
again indicating possible variability in the Fe Ka emission. 

Based on this comparison, the same absorbed power 
law with high energy cut-off model was fitted to each of the 
orbital phase resolved spectra and the spectral parameters 
measured. As with the pulse phase resolved spectral analy- 
ses, the cut-off energy and folding energy parameters were 
tied together, but with the column density and power law 
parameters free to vary independently. Additionally, an ab- 
sorption edge around ~5 keV and a Gaussian emission line 
at 6.4 keV were included to account for residuals observed in 
the low energy range of the spectra. The reduced \ 2 statistic 
obtained from this fitting process was 0.675 with 239 degrees 
of freedom. 

The neutral hydrogen column density and power law 
photon index parameters were each plotted against orbital 
phase, and compared with the PCA 2-9 keV light curve 
folded over the same period (Figure [131) . The reduced \ 2 val- 
ues obtained were 2.47 (null hypothesis probability of 4%) 
and 1.31 (null hypothesis probability of 26%) for the col- 
umn density and photon index respectively, indicating that 
both parameters are variable. The linear correlation coef- 
ficient values obtained for the column density and photon 
index compared to the folded orbital profiles were -0.78 and - 
0.84 respectively, indicating strong anti-correlations between 



both parameters and the orbital period folded flux. As with 
the pulse phase resolved spectral analyses, we were unable 
to constrain the equivalent width of the Fe Ka line suffi- 
ciently in order to determine whether it is variable over the 
orbital period. The orbital phase resolved spectra are shown 
in Figure [B3l in the Appendix. 



3.3 Super-orbital Variability 

As only 2 super-orbital cycles were sampled during the 
course of our observations, the value of the super-orbital 
period could not be accurately constrained. A s with the 
orbita l period analyses, the value obtained by IWen et al.l 
(2006) of 30.75 d was used for our analyses. The A, B and C 
(Figure [14]) PCA light curves covering all observations were 
folded over 30.75 d, with phase set at the flux minimum 
at MJD 53488. The super-orbital profile does not appear to 
vary with energy. 

Using the methods outlined above for the pulse and or- 
bital period analyses, we extracted five phase resolved spec- 
tra separated by 0.2 phase bins from the PCA data over the 
30.75 d super-orbital period. The ratios between each of the 
super-orbital phase resolved spectra and the best-fit model 
obtained from the overall spectrum were then plotted. The 
results are plotted in Figure [15] 

As with the pulse and orbital spectra, the ratio plots in- 
dicate significant spectral variability over the super-orbital 
period. Comparison with the results of the simulations pre- 
sented in Figure [8] indicates likely variability in the slope of 
the power law and in the column density parameters, with 
no evidence for variability in the high energy cut-off param- 
eters. The same features around ~5 keV and ~6 keV as 
observed in the pulse and orbital plots are present, which 
we again attribute to the xenon L-edge and fluorescent iron 
emission line. 

Based on this comparison, the same absorbed power law 
with high energy cut-off model was fitted to each of the spec- 
tra and the spectral parameters measured. Exactly the same 
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Figure 10. Variability of the column density (top) and power 
law photon index (bottom) over the pulse period, with error bars 
representing 90% confidence intervals. The histogram represents 
the RXTE PCA 2-9 keV light curve folded over the best fit 2.65 
h pulse period. The red dashed lines indicate the mean column 
density and photon index. 
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method was used as was with the pulse and orbital phase 
resolved analyses. The reduced x 2 statistic obtained from 
this fitting process was 0.6f5 with 238 degrees of freedom. 

The neutral hydrogen column density and power law 
photon index parameters were then plotted against super- 
orbital phase, and compared with the PCA 2-9 keV light 
curve folded over the same period (Figure [To) . The reduced 
X 2 values obtained were f.88 (null hypothesis probability 
of 11%) and 1.78 (null hypothesis probability of 13%) for 
the column density and photon index respectively, indicat- 
ing that both parameters are marginally variable. The linear 
correlation coefficient values obtained for the column den- 
sity and photon index compared to the folded super-orbital 
profiles were -0.40 and -0.84 respectively, indicating a weak 
and strong anti-correlation respectively between both pa- 
rameters and the super-orbital period folded flux. As with 
the pulse and orbital analyses, we were unable to constrain 
the equivalent width of the Fe Kq line sufficiently in order 
to determine whether it is variable over the super-orbital 
period. The super-orbital phase resolved spectra are shown 
in Figure [B4l in the Appendix. 



3.4 Other Variability 

A search for periodic variability between 500 - 1,500 s was 
performed in order to determine whether the flaring episodes 



Figure 11. PCA Standard2 A (top), B (middle) and C (bottom) 
light curves of all observations folded over the 11.599 d orbital 
period, with error bars determined from the standard error. 
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Figure 12. Ratios between the orbital phase resolved spectra and 
the best fit model from the overall spectrum (top inset figures) 
plotted against the folded orbital profile (histogram). 



10 S. A. Farrell et al. 



E 4 




CL 



i. 



T 



S I s 



0.5 



1.5 



Orbital Phase 



Figure 13. Variability of the column density (top) and power law 
photon index (bottom) over the orbital period, with error bars 
representing 90% confidence intervals. The histogram represents 
the RXTE PCA 2-9 keV light curve folded over the 11.599 d 
orbital period. 



which were seen to occur during the pulse maximum were 
periodic, pseudo-periodic or stochastic in nature. The PCA 
2-9 keV (A) light curves for observations 1, 10, 11, 12, 
13, 14, 15, 20 and 22 were analysed, as each had good pulse 
phase coverage and the pulse modulation could be clearly 
seen. An efsearch analysis was performed on each of the 
light curves in the 500 - 1,500 s range with 1 s resolution. 
The strongest peaks in the x 2 vs period plots were found to 
lie between 1,100 and 1,470 s in all the data sets except for 
that pertaining to observation 11, which showed a significant 
peak at 710 s. Table [3] lists the flare recurrence timescales 
and FWHM uncertainties as measured from the x 2 vs period 
plots. 

The location and width of the strongest peak in each 
data set varied considerably with time. In order to determine 
whether the flare recurrence times varied smoothly over the 
orbital or super-orbital periods, the values in Table [4] were 
plotted against orbital (Figure I17jl and super-orbital (Fig- 
ure [18]) phase, with the 710 s value excluded as an outlier. 
No correlation between the flare recurrence times and the 
orbital or super-orbital period folded flux profiles was ob- 
served, indicating that the recurrence timescales of the ob- 
served flares do not vary periodically over either the orbital 
or super-orbital period. 

It is interesting to note, however, that no flares were 
observed during the minimum flux of the super-orbital cy- 
cle. The flares were observed in conjunction with the max- 




Super-orbital Phase 

Figure 14. PCA Standard2 A (top), B (middle) and C (bottom) 
light curves of all observations folded over the 30.75 d super- 
orbital period, with error bars determined from the standard er- 
ror. 
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Figure 15. Ratios between the super-orbital phase resolved spec- 
tra and the best fit model from the overall spectrum (top inset fig- 
ures) plotted against the folded super-orbital profile (histogram). 
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Figure 16. Variability of the column density (top) and power 
law photon index (bottom) over the super-orbital period, with 
error bars representing 90% confidence intervals. The histogram 
represents the RXTE PCA 2-9 keV light curve folded over the 
best fit 30.75 d super-orbital period. 



Table 4. Recurrence timescales of the flaring episodes during 
pulse maximum, with FWHM uncertainties. 



Obs. Date 
(MJD) 



Period 

(s) 



f 


53506 


1,330 ± 60 


10 


53533 


1,300 ± 50 


11 


53535 


710 ± 20 


12 


53717 


1,400 ± 300 


13 


53721 


1,370 ± 80 


14 


53723 


1,190 ± 70 


15 


53726 


1,280 ± 60 


20 


53741 


1,470 ± 30 


22 


53747 


1,100 ± 40 



imum flux in the pulse cycle (phases ~0.2 - 0.7), and were 
completely absent during the pulse minimum. During the 
minimum of the super-orbital cycle, the pulse modulation 
is absent, and as such no flares were observed. We interpret 
this relationship as indicating that the super-orbital modu- 
lation is a result of variable accretion through some as-yet- 
unknown mechanism, as the pulse modulation would not be 
expected to be observed during phases of low accretion. 
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Figure 17. Flare recurrence times from Table|4] (squares) plotted 
against orbital phase. The histogram represents the RXTE PCA 
2 — 9 keV light curve folded over the 11.599 d orbital period. 
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Figure 18. Flare recurrence times from Table [4] (circles) plotted 
against super-orbital phase. The histogram represents the RXTE 
PCA 2 — 9 keV light curve folded over the 30.75 d super-orbital 
period. 



4 DISCUSSION 

A total of 22 individual pointed observations were performed 
with the RXTE, covering 22 pulse cycles, 6 orbital cycles 
and 2 super-orbital cycles. Significant temporal and spec- 
tral variability was observed from 2S 0114+650 over a range 
of timescales. The pulse, orbital and super-orbital modula- 
tions were all clearly present. An additional modulation over 
timescales of ~1,300 s was also present, coincident with the 
maximum of the pulse period. The overall 3-50 keV spec- 
trum was best fit with the same absorbed cut-off power law 
reported by other authors, albeit without the previously re- 
ported Fe Kq emission line at 6.4 keV or the CRSF at 22 
keV. The following sub-sections present a discussion of the 
major results obtained. 



4.1 Eclipses 

As described in § 3, the light curve of observation 22 (Fig- 
ure IA4|) exhibits an episode of very low flux which could 
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be indicative of an eclipse. Comparison with the ASM light 
curve folded over the 11.599 d orbital period shows that 
this episode of very low-flux does not occur at the expected 
orbital minimum but instead just after the maximum. In 
addition, while the pulse modulation is not present during 
observations 4 and 18, it can be seen clearly in observations 8 
and 14, all of which took place during similar orbital phases 
as observation 22. Thus, the low-flux episode in observation 
22 is most likely due to a temporary suspension of accretion, 
possibly due to density fluctuations in the stellar wind. 



4.2 Fe Kq Emission 

We did not find that the inclusion of an Fe Ka Gaussian 
emission feature improved the fit significantly in our analysis 
of the combined spectrum from all 22 observations. However, 
during our phase resolved spectral analyses over the pulse, 
orbital and super-orbital periods - where the 1% system- 
atic error was removed and the residuals due to the xenon 
L-edges were modelled using an absorption edge, so as to 
more accurately constrain the column density (see §2) - sig- 
nificant residuals around ~6.4 keV were observed in many 
of the spectra. Addition of an Fe Ka Gaussian emission line 
was thus required in order to obtain the best fit. The energy 
for this line was fixed at 6.4 keV for the phase resolved anal- 
yses, as the energy could not be well constrained. While the 
equivalent widths of these lines ranged between ~60 - 800 
eV, the estimated uncertainties were too large to determine 
whether the equivalent widths were correlated with any of 
the known peri o diciti es. 

lHall et all (2000) note that in their observations the 
equivalent width of the iron line could only be constrained 
during the low-state of the spin period, as the emission fea- 
ture was most strongly p resent during the episodes of low- 
flux. iMasetti et al.l (|2006h observed the same phenomenon, 
concluding that the Fe emission was overwhelmed by the 
continuum emission in the high state and was thus only de- 
tectable during the minimum flux periods between pulses. 
In contrast, the relatively poor spectral resolution and sen- 
sitivity of the RXTE PCA meant that the parameters of the 
Fe line could not be adequately constrained in our data. Fur- 
ther observations with an instrument with greater spectral 
resolution and sensitivity than the RXTE PCA are required 
before a firm conclusion can be drawn as to the nature and 
variability of the fluorescent iron emission in this system. 
Understanding the origin of this emission could lead to a 
clearer picture of the geometry and dynamics of the system, 
potentially supplying further evidence for the presence of a 
super-slow pulsar and providing an answer to the question 
of the origin of the super-orbital modulation. 



4.3 Cyclotron Resonance Scattering Features 

The C RSF at ~22 keV reported by iBonning fe Falangal 
(2005) w as not evident in any of ou r spectra, support i ng the 
results of lMasetti et al.l l|2006h and lden Hartog et all (|2006l ) 
who also failed to detect this feature. 

If the magnetic field of th e neutron star in 2S 01 1 4+65 
is > 10 12 G as suggested bv iLi fe van den Huevell ljl999h . 
we would expect a CRSF in the spectrum at > 9 keV. 
iHeindl et al.l (|2004h showed there to be a rough correlation 



between the energy of the cut-off and the cyclotron line en- 
ergy in the spectra of accreting X-ray pulsars (see their Fig- 
ure 4). Using a value of ^6 keV for the cut-off energy as 
determined from our spectral fits, 2S 0114+650 should have 
a CRSF line energy of ~6 keV and therefore a magnetic field 
strength of ~10 12 G. 

The scattering cross section has a strong dependence 
on the angle to the magnetic field, so that CRSFs provide 
an excellent diagnostic of the emission geo metry and the 
phys ical conditions of the radiating plasma (jHarding fe Lail 
120061 ) . The angle between the photon direction and the mag- 
netic field vector is thus a determining facto r in the detection 
of cy clotron features in the spectrum (e.g. ISchonherr et al.l 
2007). The shape and relative depths of the CRSF first and 
higher harmonics are thus highly dependent on the view- 
ing geometry and th e geometry of the scattering plasma 
IjHarding fe Laill2006T ). As a result, variations in the features 
are expected with the varying viewing geometry produced 
by the rotation of the neutron star. For the case where the 
photon direction and magnetic field vector are aligned, no 
cyclotron resonance scattering would be expected. In this 
scenario it may be possible to place some constraints on 
the accretion geometry (e.g. fan beam/cyclinder vs pen- 
cil beam/slab geometry), rotation patterns and accretion 
rate through the abs ence of observable CRSF lines (e.g. 
ISchonherr et al.ll2007f l. 

Alternatively, the absence of any CRSF features in 
the X-ray spectrum could indicate that the magnetic field 
strength is outside the observable bandwidth, indicating a 
field strength of either < 10 12 G or > 10 12 G. In either of 
these scenarios the CRSF would occur at energies too low 
or too high for us to observe in our data. 



4.4 The High-Energy Tail 

The detection of a 'high energy tail' in the spectrum of 2S 
0114+650 from analysis of INTEGRAL IBIS-ISGRI data 
was reported bv lden Hartog et al.l (|2006l ). which required a 
combination of thermal bremsstrahlung and power law mod- 
els to describe the spectrum at high en ergies. The i r find - 
ings are thus i ncons is tent with those of Hall et alj |2000), 
IMasetti et ail (120061 ). iMukheriee fe Paull (|2006h and with 
our own results. 

In order to test the validity of their model, we at- 
tempted to fit a thermal bremsstrahlung + power law model 
to our combined PCA + HEXTE spectra in the 3-50 
keV range. While this model provided an acceptable fit, 
significant residuals were seen at ~5 keV indicating that 
the model does not accurately describe the spectrum. Due 
to the poor statistics in the HEXTE data above 50 keV 
we could not confirm the existence of a tail in the spec- 
trum > 70 keV. However, the poor residuals at low energies 
obtained with the combined thermal bremsstrahlung/power 
law model raises questions about the presence of a high en- 
ergy tail i n the spectrum fr o m thi s system. It should be 
noted that Iden Hartog et aLI |2006l ) also found an accept- 
able fit to the spectrum using a single power law model 
without the inclusion of a high energy exponential cut-off. 

Additionally, in other binary systems the appearance of 
a high energy tail has been linked to the onset of compact 
jets. Radio observations of 2S 0114+650 performed in the 
240 - 1400 MHz range with the Giant Metrewave Radio Tele- 
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scope found no evidence for any radio emissio n, thus argu- 
ing a gainst the presence of jets in this system (jFarrell et al.l 
12007t h 



4.5 The ~2.7 hr Pulse Period 

From our analysis of the PCA timing data, we have ob- 
tained a new value for the pulse period at 2.65 ± 0.003 
h. This value is inconsis tent with the value reported by 
iBonning fc Falangal (|2005l ). indicating that the neutron star 
is continuing to spin-up over time (see lHall et all |2000| ; 
ISood et al.ll2006l ). Separate analysis of the two observation 
runs showed no significant difference between the values of 
the pulse period during the six months separating the runs, 
providing an upper limit for P spin I P spin during MJD 53505 
- 53747 of < - 4.6 x 10~ 3 yr _1 . Our results are thus con- 



Our analysis of the phase resolved spectra showed 
marginal evidence for variability in the power law photon 
index over the pulse period, anti-correlated with the flux. A 
similar anti-correlation between the photon in dex and the 
pulse flux has been observed in Her X-l (e.g. ISoong et al.l 
Il990l ). The observed similarities between the spectral vari- 
ability of 2S 0114+650 and Her X-l over the pulse periods 
could thus strengthen the case that the ~2.7 h period is 
indicative of the neutron star spin period. However, the un- 
certainties for the photon index values are too large to draw 
a firm conclusion either way. No significant correlated vari- 
ability of the neutral hydrogen column density over the pulse 
period could be determined. 

Usin g the derived values of the w ind velocity (1,200 



Corral fc Koenigsbergerl Il987). orbital tangential 



km s 

velocity (~271 km s 



assuming an almost-circular orbit), 



sistent with the value of 
lHall et al.l <|2000h . 

Our results also indic ate that the model proposed by 
iKoenigsberger et al.l (|2006l ) - which ascribed the ~2 h pulse 
period to tidally-induced pulsations of the B-supergiant 
companion - is inconsistent with the observed properties 
of this source. They showed that tidal oscillations of the 
companion star's surface could occur at both the ~2 h and 
~30 d timescales. However, their model predicts that the 
~2.7 h period should vary significantly on short times-scales 
between ~1.6 - 3.4 h (depending on the system parame- 
ters). In addition, the link between tidal oscillations of the 
donor star and periodic ma ss ejection episodes is still unclear 
|Koenigsberger et al.ll2006h . Our observations thus indicate 
that this model does not adequately explain the mechanism 
behind the pulse period, as the stability of this period has 
been demonstrated over timescales of days to months. 

The pulse profile shown in Figure [7] is approximately 
sinusoidal with features around phase 0.2 and 0.85, with 
the latter possibly arising fro m the edge of the em ission 
region at the second pole fe.g. IWang fc Welter! Il981). The 
profile of the modulation does not appear to vary signifi- 
cantly with energy, implying an absence of cyclotron scat- 
tering (as any cyclotron features in the spectrum would be 
expected to vary wit h the rotation of the magnetic field; e.g. 
lHarding fc La] |2006). supporting the results of our spectral 
fitting. In addition, while no clear-cut dependence of the 
pulse profile on the pulse p eriod is apparent for X-ray pul- 
sars, IWang fc Welterl l|l98ll ) noted that many slow pulsars 
{Pepin > 100 s) show relatively smooth and sinusoidal pro- 
files. Our results thus continue to support a pulsar nature 
for the compact object in 2S 0114+650. 

The approximately sinusoidal profile and the lack 
of change with energy is consistent with the profiles of 
other X-ray pulsars with lum inosities below ~10 36 erg s _1 
l| White. Swank, fc Holtl Il983h . Previous attempts to con- 
strain the distance to 2S 0114+650 have yielded a range of 
value s (1.8 ± 0.2 kpc, ~2-5 kp c, ~3.8 kpc, ^3 kpc; Aab et ahl 
19831: ICrampton et all Il985l: IChevalier fc Ilovaiskvl 1 19981 : 



2 x 10 -3 yr _1 obtained by orbital angular velocity (6.3 x 10 rad s ), and typical 



Koenigsberger et al Kooa res pectiv ely) , with the most thor- 
ough attempt bv lReig et al.l l|l996h providing a poorly con- 
strained distance of 7.0 ± 3.6 kpc. Taking 10 36 erg s _1 as 
the luminosity up per-limit in the 0.5 - 60 keV band from 
IWhite et al.l (|l983l ) we can provide further constr aints on the 
distan ce. Assuming a lower limit of 3.4 kpc from lReig et al.l 
(1996), we derive a new distance estimate of 4.5 ± 1.5 kpc. 



values for the other parameters l|lkhsanovll2007h . we calcu- 
lated the M in order to achieve the observed average spin-up 
rate. The results of these calculations yielded values of ~9 
x 10" 9 M Q yr" 1 for a neutron star and ~2 x 10" 3 M Q 
yr _1 for a white dwarf. The derived mass accretion rate for 
a neutron star is approximately two orders of magnitude 
les s than the mass loss ra te for Bl la supergiants derived 
bv lBarlow fc Cohenl d 19771) from infrared observations and, 
more recently, by Crowther. Lennon fc Walbornl l|2006l ) from 
optical observations. 

The mass accretion rate would be expected to be much 
less than the mass loss rate of the donor star due to ineffi- 
cient capture of material from the stellar wind. Our derived 
mass accretion rate is thus consistent with the expected mass 
loss rate of the supergiant. In contrast, the derived mass 
accretion rate for a white dwarf is three orders of magni- 
tude greater than the expected supergiant mass loss rate, 
clearly excluding a white dwarf as the compact object in 
this system. In light o f this, the evolution of the pulse pe- 
riod ((Sood et al.ll2006h . and the spectral similarities shared 
with other binary X-ray pulsars, we conclude that the ~2.7 
h period must indeed be the spin-period of a neutron star. 

The recent rep orts by iMattana et al.l l|2006l ) and 
|Pe Luca et al.l (120061 ) of two other candidate pulsars with 
hour-long spin-periods (~5 h for 4U 1954+319 and ~6.7 h 
for IE 161348-5055.1 respectively) - if confirmed - would ap- 
pear to indicate that super-slow rotating neutron stars are 
more prevalent than originally suspected. Ilkhsanovl (120071 ) 
has shown that the assumption of a super-critical initial 
magnetic field strength (i.e. B 3> 10 13 G) is not necessary 
if the propeller phase consisted of both supersonic and sub- 
sonic propeller phases. Rapid spin-down would occur over 
very short timescales due to interactions between the ac- 
creting material and the neutron star magnetic field. In this 
setting, even spherical accretion from a stellar wind would 
result in significant angular moment um transfer due to the 
orbital motion of the neutron star ijlkhsanovl 120071 '). lead- 
ing to spin-down timescales consistent with the projected 
lifet ime of the donor star. 

Ilkh" sanov (I2007D showed that in the Davies-Fabian- 
Pringle l|Davies. Fabian, fc Pringie! 1 19791) scenario a mass 
accretion rate of ~10 14 g s _1 is required to spin 2S 0114+650 
down to a period of ^10 4 s. This derived accretion rate is a 
factor of 30 smaller than the current mass transfer rate in- 
ferred from the X-ray luminosity of the source. However, the 
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evolutionary tracks of massive stars predict a main sequence 
09.5 V progenitor star for a Bl la supergiant in a previous 
epoch ()Mevnet et al.l[l994h . with a mass outflow rate of 30 - 
40 times smaller than the donor star in 2S 0114+650 inferred 
from optical observations. Hence, the mass accretion rate in 
a previous (spin-down) epoch would be ~10 14 g s^ 1 , consis- 
tent with the value required by the Davies-Fabian-Pringle 

scenario. 

The pr oblem common t o both the lLi fc van den HueveJ 

l|l999l ) and llkhsanovl (|2007r i scenarios is that the equilib- 
rium period for 2S 0114+650 is predicted to be < 26 min, 
which the spin-period will rapidly approach on timescales 
of < 100 yr fo r disc accretion and < 1,000 yr for stellar 
wind accretion l|lkhsanovll2007l) . The probability of observ- 
ing the pulsar in 2S 0114+650 during the long-period phase 
is thus almost negligible, indicating either 2S 0114+650 is a 
young accretor (and we have beaten all the odds to observe 
it in this stage) or our current understanding of pulsar evo- 
lution is deficient. If the pulse period continues t o decrease 
at th e current average rate of ~ —2 x 10~ 3 yr -1 (Ha ll et all 
120001 ). the pulsar will reach its equilibrium period in < 400 
yr. This timescale is thus consistent with an infant X-ray 
pulsar spinning-up through stellar wind accretion. 




ORBITAL PHASE 



Figure 19. Variability of the photon index over orbital phase 
of EXO 1722-363. Phase 0.5 corresponds to the maximum of the 
orbital period, with phases 0.9 - 1.1 correspondin g to the phase 
of eclipse (data taken from IThompson et alj|2007t) . 



ties in the orbital properties. We thus ascribe the observed 
spectral variability to 2S 0114+650 passing behind a heavily 
absorbing region close to the base of the wind of the donor 
star. 



4.6 The ~11.6 d Orbital Period 

The observed variability of the neutral hydrogen column 
density over the orbital period indicates that variable ab- 
sorption by the dense stellar wind throughout the orbit is 
the mechanism behin d the orbital modulation , as su ggested 
by lHall et all ||2000D and iGrundstrom et all (|2007h . Simi- 
lar variability in the column density has been obs erved in 
the HMXB pulsars Vela X- l and EXO 1722-363 (jNaeasd 
ll99ll : lThompson et alJl2007t ). The latter system in particu- 
lar provides a good comparison, as both the column density 
and the power law photon index were observed to vary over 
the orbital period, anti-correlated with flux (see Figure [19] 
for a plot of the photon index over orb ital phase from the 
data presented in lThompson et al.ll2007r i . EXO 1722-363 is a 
highly obscured binary system that exhibits regular eclipses. 
The observed increase in the column density and photon in- 
dex during eclipse are thus a result of absorption effects. 

It is not possible to determine from our results whether 
accretion in an unperturbed stellar wind is sufficient to ac- 
count for the observed changes in the column density. Thus, 
we cannot determine whether structures such as accretion 
wak es (similar to those detected in Vela X-l and EX O 1722- 
363; lHaberl fc White! Il990l ; IThompson et al.l 120071 ') are re- 
quired to account for the observed absorption. However, the 
shape of the column density variations with orbital phase is 
consistent with an undisturbed wind, and not like the highly 
disturbed win d models obtained from simulations (see Fig- 
ures 2 and 8 of lBlondin. Stevens, fc Kallmanlll99 if ). Further 
observations and modelling are required before the geometry 
of the wind can be constrained. 

While our results do not exclude periodic eclipses of 
the neutron star by the supergiant companion, we found no 
strong evidence for eclipses in the PCA light curves. We thus 
conclude that 2S 0114+650 does not exhibit t he properties 
of an eclipsing binary system as suggested bv lCorbet et al.l 
(1999). However, the similarities between the spectral vari- 
ability in 2S 0114+650 and EXO 1722-363 hint at similari- 



4.7 The ~30.7 d Super-orbital Period 

Her X-l is the prototypical system in which the warped, pro- 
cessing disc scenario has been shown to explain the observed 
super-orbital modulation very well. The column density in 
this system has been observed to vary out of phase with the 
super-orbital cycle, due to variable obs curation of the cen - 
tral source by the precessing disc warp (Na ik fc Pauli r2003). 
The lack of variability of the column density over the super- 
orbital period in 2S 0114+650 thus indicates that this mod- 
ulation is not due to local changes in the column density, 
and therefore precludes variable absorption by a warped, 
precessing accretion disc as the cause of the super-orbital 
modulation. 

The lack of significant Fe Ka emission precludes vari- 
able reflection by a precessing disc, as strong iron emis- 
sion due to reprocessing of X-rays in a disc would be ex- 
pected to be present and variable in this scenario (e.g. 
iGeorge fc Fabianlll99ll ). In addition, the apparent increase 
in the power law photon index during the super-orbital min- 
imum is in exact opposition to the behaviour of two of the 
othe r well studied super -orbital systems - LMC X-4 and Her 
X-l (|Naik fc Paulll2003l ') - both of which exhibit photon in- 
dex variability modulated over the super-orbital period. In 
both of these systems the super-orbital modulation is at- 
tributed to the warped precessing disc model, once again 
showing that the 30.75 d period in 2S 0114+650 does not fit 
the standard model. 

If the super-orbital modulation was tied to accretion 
onto another compact object (i.e. a triple system), we would 
expect to see variability of the column density and possibly 
the photon index in a similar fashion to the orbital mecha- 
nism discussed above. While modulation of the eccentricity 
of the 11.6 d orbit by a third massive body would not lead 
to variability in the column density, this model would yield 
a long-period of S> 30.75 d. 

Similarly, if precession of the neutron star-spin axis were 
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the cause of the super-orbital modulation, the precession 
period would be inversely proportional to the spin period 
i|Wassermanll2003h . We would thus expect to see the super- 
orbital period vary, anti-correlated with the spin-up of the 
neutron star . The apparent lon g-term stability of this period 
as shown bv lSood et al.l (|2006h indicates that this latter sce- 
nario is unlikely. In addition, the expected precession period 
for a neutron star spinning with a ~2.7 hr period is 3> 30.75 
d. 

The periodic episodes of enhanced mass accre- 
tion in the tidally-driven oscillation model proposed by 
iKoenigsberger et alj l|2006h would n o t be e xpected to cause 
column density variability. iHuded (|l978r ) analysed 1,250 
photographic plates taken of 2S 0114+650 between 1928 and 
1977 at optical wavelengths and could not discern any sig- 
nificant variability in the optical magnitude greater than 0.3 
mag over any timescale. Pulsations in the supergiant donor 
star should result in significant variable optical luminosity 
thus making this scenario unlikely. However, it should be 
noted that no long-term optical monitoring has been per- 
formed since 1999 when the super-orbital period was first 
detectable in the ASM data, preventing us from ruling out 
this model as an explanation for the super-orbital variability. 

Another possibility is that the neutron star under- 
goes periods of enhanced accretion over the super-orbital 
timescale due to the presence of another donor star in an ec- 
centric orbit of ~30 d. As the second donor star approaches 
periastron passage, accretion onto the neutron star would 
commence through Roche lobe overflow. The increase in flux 
over the super-orbital period would thus be due to phase 
dependent episodes of Roche lobe overflow, which would 
be associated with the formation of a transient disc dur- 
ing very specific super-orbital phases. This would result in 
an increase in the accretion rate and therefore an increase 
in flux, but would not necessarily result in variability of the 
observed spectral parameters. Previous optical observations 
have identified only one candidate for the optical counter- 
part within the X-ray error circle. If the super-orbital mod- 
ulation is due to the presence of a second donor star in a 
triple system, the third star would have to be low-mass with 
very low luminosity. 



4.8 Other Variability 

Pseudo-periodic flaring episodes were observed around 
~1,300 s at the maximum pulse phase, with the recurrence 
timescale of these flares seen to vary significantly over 
time. No evidence for a correlation between the flare 
recurrence timescales and either the orbital or super-orbital 
period flux was found. When taken into context with the 
894, 850 and ^2,0 s p er iods previously r eporte d by 
Kocnigsberger et all (|l983h . lYamauchi et~al] (|l990h and 
Apparao et al.l l|l99ll) respe c tively and the X-ray flickering 
reported by Masetti et al.l l|2006 ), we conclude that the 
observed ~1,300 s flares are most likely due to an aperiodic 
process with a stochastic recurrence timescale. Possible 
explanations include "clumping" of the material in the 
accretion stream due to turbulent processes, or possibly 
feedback between the radiation pressure from the emission 
region and the accreting material. Both explanations imply 
a correlation between the flare recurrence timescales and 



the mass accretion rate. 



5 CONCLUSIONS 



We have performed an extensive study of the X-ray vari- 
ability of 2S 0114+650 with the RXTE satellite. The good 
spectral coverage provided by the PCA and HEXTE instru- 
ments have allowed us to undertake the first in-depth wide- 
band study of the super-orbital modulation. The overall X- 
ray spectrum in the 3-50 keV range was found to be well 
described by the standard absorbed power law model with 
a high energy exponential cut-off, with spectral parameters 
consistent with those reported by other authors from other 
instruments and epochs. 

Marginal evidence for variability in the power law pho- 
ton index over the pulse period was identified, which could 
be attributed to the varying viewing geometry of the accre- 
tion region with the spin of the neutron star. Similar vari- 
ability has been observed from other X-ray pulsars such as 
Her X-l, supporting the conclusion that 2S 0114+650 con- 
tains a super-slow X-ray pulsar. The near-sinusoidal pulse 
profile and the lack of change with energy puts 2S 0114+650 
in the class of low luminosity (< 10 36 erg s _1 ) X-ray pulsars, 
allowing us to provide a new distance estimate of 4.5 ± 1.5 
kpc. 

Variability of the neutral hydrogen column density over 
the orbital period was observed, indicating that the line- 
of-sight motion of the neutron star through the dense cir- 
cumstellar environment is the mechanism behind the orbital 
modulation. Additional variability in the power law photon 
index anti-correlated with flux over the orbital period was 
identified, which we speculate is due to absorption effects as 
the neutron star passes behind a heavily absorbing region 
near the base of the supergiant companion's wind. 

In contrast, no significant variability of the column den- 
sity was observed over the super-orbital period. While super- 
orbital periods are commonly attributed to variable obscura- 
tion by a precessing warp in an accretion disc, the absence of 
super-orbital column density variability indicates that this 
is not the mechanism behind the super-orbital modulation. 
Similarly, the absence of significant iron emission or reflec- 
tion components in the spectrum argues against variable 
reflection by a precessing warped disc. While the slope of 
the spectrum was relatively stable throughout the super- 
orbital cycle, a significant increase in the photon index was 
observed during the minimum phase. We conclude that the 
super-orbital period is tied to variability in the mass accre- 
tion rate due to some as yet unidentified mechanism. 
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APPENDIX A: INDIVIDUAL LIGHT CURVES 

The individual light curves for the 2-9 keV (A) energy 
band from each of the 22 observations are included in this 
Appendix. 



APPENDIX B: PHASE RESOLVED SPECTRA 

The spectra from the pulse, orbital and super-orbital phase 
resolved analyses are included in this Appendix. 
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Figure Al. The 2-9 keV PCA light curves for observations 1 to 6, with the time axes showing hours elapsed since MJD 53505. 
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Figure A2. The 2-9 keV PCA light curves for observations 7 to 11, with the time axes showing hours elapsed since MJD 53505. 
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Figure A3. The 2-9 keV PCA light curves for observations 12 to 17, with the time axes showing hours elapsed since MJD 53715. 
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Figure A4. The 2-9 keV PCA light curves for observations 17 to 22, with the time axes showing hours elapsed since MJD 53715. 
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Figure B4. PCA 3-30 keV spectra from the super-orbital phase resolved spectral analyses ordered from bottom to top at 5 keV: 0, 
0.8, 0.2, 0.6, 0.4. 



